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EMPIRICAL FORMULA FOR THE VOLUME OF 
ATMOSPHERIC AIR. 



BY a. W. HILL. 

The formula of Mariotte and Gay-Lussac is generally employed, in the 
laboratory, to reduce volumes, observed under one tension and temperature, 
to those which would have place under other tensions and temperatures. 
But Regnault, about 1845, made several series of experiments, which, if 
they may be relied upon, establish marked deviations from this formula. 
These experiments are detailed in the Me moires de I'Acade'mie des Sciences 
de Paris, Tom. XXI. I propose to investigate a modification of the for- 
mula, the introduction of which makes it possible to satisfy nearly these 
experiments. 

Let T denote the temperature, here always expressed in degrees of the 
centigrade scale; P the tension or pressure, measured by the altitude, in 
metres, of a column of mercury, it is capable of supporting; the mercury 
being at the temperature and under the action of gravity which obtains 
at Eegnault's laboratory; and let V denote the volume. Then, for any 
given mass of air, these three quantities are so connected that if any two of 
them are assigned the remaining third is immediately determined. That 
is, we must have function ( V, P, T) = 0, 

or, solved with respect to V, 

V= function (P, T). 
But the mode, in which T is to be measured, is arbitrary, and we may take 
atmospheric air as the thermometric substance, and assume that T increases, 
in direct proportion, as the volume, under constant pressure, increases. 
This gives 

7= F(P) + f(P).T. 
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It is here taken for granted that, whatever may be the density of the air 
inclosed in the thermometer, its indications will be the same. It is true 
that the usual custom of experimenters has been to measure temperatures 
by the augmentation of tensions under constant volume; but, when Mari- 
otte's law holds, this gives results identical with those obtained by the for- 
mer method. In this case we shall have to write the equation 

P = F(V) + f(V).T. 
But the first equation seems preferable. 

Now since, for any given constant temperature, the volume ought to be a 
function of the tension similar to what it is at any other temperature, it fol- 
lows that, if F[P) is supposed to consist of a series of terms, each of the 
form KP k , where K and h are constants, so that we may write 

F{P) = E.KP", 
then we ought to have 

f(P) = I.K X P", 
where K x denotes a constant, in general, different from K. Thus we should 
have 

V=2.\K+K 1 T~\P k . 

The formula of Mariotte and Gay-Lussac assumes that I\P) and /(P) 
contain each only one term, in which h = — 1. But Regnault's experi- 
ments having shown the insufficiency of this, it is in order to see whether 
agreement between theory and observation cannot be brought about by an- 
nexing to V an additional term, for which k has a value different from — 1. 
Thus let us suppose that 

V= [K+K^P- 1 + [K'+K^TIPP- 1 

_ K+ K'PP , K X +K{P P T 
p + p J- 

As V contains a factor, which is directly proportional to the mass of air 
considered and inversely as the unit assumed for the measurement of vol- 
umes, we prefer to write the preceding equations thus 

V=K[l±^ + a -±£Tpq 
When the temperature is constant, the volumes are represented by the formula 
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that is, the resultfroin Mariotte's law must be multiplied bythe factor 1-f aPP, 
which differs but little from unity. « is a small constant which measures 
the amplitude of the deviations from this law ; while /? is a constant expo- 
nent, which is so chosen that the more or less rapid variation of the devia- 
tions, in passing from one tension to another, may be represented as well as 
possible. It is evident that, in this manner, we get the utmost advantage 
that can be derived from the addition of a single term to V. 

The experiments of Eegnault may be divided into two classes; first, those 
where, the temperature remaining nearly constant, the volumes of the same 
mass of air, under different pressures, were observed; second, those where, 
the volumes remaining nearly the same, the tensions were observed at the 
temperature of freezing and boiling water. It is obvious that experiments 
of these two kinds extended over a sufficient range of tension would afford 
the data requisite for obtaining the values of the four constants a, a', a" 
and /3 which enter into our adopted formula. 

The experiments of the first class are enumerated at pp. 374 — 379 of the 
volume quoted above. As the temperature is nearly the same for all, we 
assume that they have been made at the average of all the noted tempera- 
tures which is 4 .747. 

To save labor we may take the average of the observed volumes and ten- 
sions when they are nearly alike. In this way Regnault's 66 experiments 
are reduced to the 23 given in the following table. It may be noted that 
Fis here expressed by the number of grammes of mercury required to fill 
the volume. The column containing log (P V) exhibits the deviation from 
Mariotte's law; did this law exactly hold, the numbers in this column would 
be identical for each series. It will be noted that, in general, they diminish 
with increasing pressures. The volumes being supposed to be represented 
by the equation 

a preliminary investigation has given the approximate values 
a = — 0.0024337, £ == 0.645. 

With these have been computed the values of the expressions which stand 
at the head of the two last columns of the table, and which serve to obtain 
the coefficients of the equations of condition to be given presently. 

As the mass of air operated on was different in each series of experiments, 
K will have 9 different values; it can, however, be eliminated. Taking the 
common logarithms of each member of the equation last given, 
logiJT -f log(l + aP) =P log(PF). 
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First Series -j 

Sec'd Series ■< 

Third Series 
F'rth Seri 

Fifth Series < 
Sixth Series - 
SVth Series « 



Eig'h Series < 
Ninth Series < 



V. 

1939.76 
969.65 

1939.37 
970.40 

642.82 

1939.72 
969.78 

1940.65 
979.78 

1939.85 
970.29 
626.91 

"1940.23 
970.32 
685.11 
675.15 

' 1941.23 
969.98 
633.82 

1940.44 
970.53 

1945.06 

1053.78 



i>.N„.o.».i„ g (pn^. I |^iogp. 

0. m 73899 4 3.156387 0.8244 —0.1083 

1. 47630 4 3.155790 1.2897 +0.2182 

2. 11228 3 3.612412 1.6262 0.5281 
4. 21020 3 3.611254 2.5430 1.5876 
6. 35034 2 3.610886 3.3213 2.6664 

2. 06887 3 3.603472 1.6045 0.5066 

4. 12663 6 3.602268 2.5102 1.5452 

4. 14235 2 3.905194 2.5164 1.5532 

8. 17850 3 3.903803 3.9155 3.5737 

4. 21910 4 3.912988 2.5465 1.5921 

8. 40648 4 3.911516 3.9863 3.6857 

12. 98195 1 3.910545 5.2926 5.8925 

6. 77001 3 4.118444 3.4623 2.8758 

13. 47353 4 4.116396 5.4226 6.1247 
19. 00213 1 4.114562 6.7913 8.6846 
19. 30191 2 4.115000 6.8612 8.8206 

6. 39003 2 4.093580 3.3347 2.6861 

12. 72859 2 4.091543 5.2248 5.7721 

19. 39954 1 4.089757 6.8842 8.8654 

9. 33401 3 4.257968 4.2676 4.1398 
18. 54702 5 4.255283 6.6842 8.4774 

47357 2 4.348632 4.8824 5.1740 

05700 2 4.346146 7.2643 9.6137 



11. 
21. 



To reduce the matter within the treatment of the method of least squares, 
it will be necessary to make some assumption regarding the probable errors 
of the observed P and V. We will, for convenience, suppose that they are 
such that the function log (P V) has a probable error equal for all the obser- 
vations; an assumption somewhat precarious, it is true, but it seems that 
we cannot easily do better. 

Let the small corrections, which it is necessary to apply to the approxi- 
mate values of log K, a. and /?, be denoted by d log K, da and <J/3, and let us 
put d log K = x, Mda = y, adfi — z, 

where M denotes the modulus of common logarithms. <Hog(PF) being 
the excess of observed over calculated log (P V), we shall have the equation 
of condition 

p/3 pp 

x + _£-_«+— ^-- R \ogP.z = Slog(PV). 



1+aPP 



1+aPP' 
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A little consideration will show that x will be eliminated by taking the 
difference of every two equations of condition arising from the same series, 
and attributing the weight ww' -5- I . w to the resulting equation, w and w' 
denoting the weights of the equations whose difference is taken, and 2" . w 
the sum of the weights of all the equations in the series. Since the coeffi- 
cients of y and z, in the equations, are all positive and nearly proportional, 
it will be advantageous to adopt a new unknown u, such that 

y — u — |-z. 

Then the equations, with the weights that ought to be attributed them, are 

1st Series 0.4653m — 0.4490z = —0.000106 Weight, 2 

(0.9168 —0.4685 == —0.000193 " £ 

2nd " 1 1.6951 —0.6869 = +0.000255 f 

(0.7783 —0.2184 = -fO.000448 § 

3rd " 0.9057 —0.4709 = —0.000252 2 

4th " 1.3991 —0.3113 = +0.000076 1.2 

f 1.4398 —0.3061 = +0.000038 ^ 

5th " 1 2.7461 —0.2764 = +0.000432 f 

(1.3063 +0.0296 = +0.000394 | 

f 1.9603 —0.0183 = +0.000002 1.2 

| 3.3290 +0.2605 = —0.000407 0.3 

»., « J 3.3989 +0.2800 = +0.000104 0.6 

om K 1.3687 +0.2787 = —0.000409 0.4 

1.4386 +0.2982 = +0.000102 0.8 

0.0699 +0.0195 = +0.000511 0.2 

(1.8901 —0.0642 = —0.000059 0.8 

7th " 1 3.5495 +0.2635 = —0.000117 0.4 

(1.6594 +0.3276 = —0.000058 0.4 

8th " 2.4166 +0.3099 = —0.000164 J^ 

9th " 2.3819 +0.4699 = —0.000005 1 

The derived normal equations are 

58.672w —0.07902 = —0.0005252, 
— 0.079m +2.44532 s= —0.0000157. 
Whence 

u = —0.000008962, z = —0.000006707, y = +0.000002216. 
da = +0.0000051, dp = +0.00276. 
Applying these corrections to the approximate values of a and /?, we get 
the final values 

a = —0.0024286, j8 = +0.64776. 

How well the experiments are represented by the formula, with these 
values of the constants, will best be seen from the following comparison of 
the values of T / P -4-F I P 1 given by Regnault and those computed from 
the formula. 
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Obs. Cal. Diff. Obs. Cal. Diff. 

1.001414 1.001133 +281 1.005437 1.006694 —1257 
1.001448 1.001132 +316 1.005703 1.006694 — 991 

iZw iZlll t 28 9 8 L004286 L004777 ~ 491 

1.001421 1.001133 +288 1.004512 1.004770 — 258 

1.002765 1.002233 +532 1.004599 1.004779 — 180 

1.002759 1.002234 +525 1.004580 1.004771 — 191 

1.002503 1.002236 +267 1.008536 1.008106 + 430 

1.003539 1.004134 —595 1.008813 1.008108 + 705 

1.003452 1.004133 —681 1.008016 1.008286 — 270 

1.003309 1.004133 —824 1.008064 1.008269 — 205 

1.002709 1.002209 +500 1-007980 1.008288 — 308 

1.002724 1.002207 +517 1.004611 1.004601 + 10 

1.002713 1.002206 +507 1.004752 1.004601 + 151 

1.002528 1.002211 +317 1.008930 1.008648 + 282 

1.002898 1.002203 +695 1.008755 1.008642 + 113 

1.002762 1.002203 +559 L006366 L005876 + 49Q 

1.003253 1.003417 —164 1.006132 1.005880 + 252 

1.003090 1.003411 —321 1.006010 1.005869 + 141 

1.003302 1.003407 —105 1.006346 1.005878 + 468 

1.003336 1.003506 —170 1.005619 1.005738 — 121 

1.003495 1.003508 — 13 1.005622 1.005736 — 114 

1.003335 1.003508 —173 1.005902 1.005832 + 70 

1.003448 1.003509 — 61 

It will be seen that the differences, in the extreme cases, amount to a fourth 
part of the observed deviation from the laAv of Mariotte. Morover it is plain 
that some cause, which varied from series to series, has operated to vitiate 
these experiments, since it is possible to determine a and /3 so that any two 
series are well represented, but not possible when all the series are included 
in the investigation. It may be noted also that the experiments in which 
the original volume was reduced to one third are not in general concordant 
with those where the reduction was to one half. 

That these discrepancies are unavoidable will be evident from the follow- 
ing exposition. Let us put 

com. log (PV) = F{P). 
The observations of Regnault may be resume'ed in the following 9 results, 
all formed by combining tolerably concordant data, 

1 2^(1.476)— 2\0.739) = 0.000598, 

2 2^4.168)— 2^2.091) = 0.001181, 

3 2(6.350)— 2 , (2.112) == 0.001526, 

4 ..... . 2 , (8.292)— 2 7 (4.182) = 0.001437, 
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5 JF(12.982)— F{ 4.219) = 0.002443, 

6 2?(13.101)— 2\ 6.580) = 0.002042, 

7 JF(19.276)- -F( 6.580) = 0.003743, 

8 F{18M7)—F( 9.334) = 0.002685, 

9 .F(21.057)— 2^(11.474) = 0.002486. 

These are the data actually furnished by Regnault for the determination of 
the function F(P). Employing the graphical method, we endeavor to con- 
struct the curve whose equation is y = F[x). One of the special values of 
F(x), as 2^(0.739), may be taken arbitrarily, and then the value of 2^(1.476) 
becomes known. This premised, we see that each of the 9 equations fur- 
nishes the length, direction and abscissae of the extremities of a chord, of the 
sought curve. Placing the chord, corresponding to the first equation, 
arbitrarily, and drawing the others on any part of the paper, but with the 
correct direction, and correct abscissae of their extremities, we endeavor, by 
imparting a motion to all their points parallel to the axis of y, to make 
them fall into line as the chords of the same continuous curve. We find 
that if we take 1, 2, 4, 6 and 7, they can be made to indicate a tolerably 
continuous curve; but then 3, 5, 8 and 9 are not satisfied. 

Again, from this graphical process, we see that there cannot be much va- 
riation of curvature between the extremities of each chord, and hence the 
tangent to the curve, corresponding to the abscissa, which is the mean of 
the abscissae of the extremities, ought to be, very approximately, parallel to 
the chord; or in other terms 

dx \ 2 / x x — x 

This gives the following values of -f-, 

dy 

x. -?-. 

dx 

1 1.108 +0.0008113 

2 3.130 0.0005686 

3 4.231 0.0003770 

4 6.237 0.0003497 

5 8.600 0.0002788 

6 9.840 0.0003131 

7 12.428 0.0002948 

8 13.940 0.0002914 

9 16.265 +0.0002594 
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From the general course of these values of dy-t-dx it may be gathered 
that this function, at first, diminishes rapidly, afterwards more slowly, and 
then tends, with higher values of x, to become nearly constant. But while 
this is the conclusion from the tout ensemble, a comparison of some of the 
values contradicts it. Thus from 1, 2 and 3, while dy-i-dx diminishes 
0.0002427 in an interval 2.0 in x, it afterwards diminishes 0.0001916 in an 
interval 1.1 of a;. All attempts then to represent these data by a curve, 
without singular points, must, evidently, show large errors. 

For the discussion of the second class of experiments, let us assume that 
a has the signification we have given it in the general formula for V. Then 
the volume remaining the same, if P and P x denote the tensions observed 
respectively at 0° and 100°, we have 

P x _ 1 + 100a' + (a+lOOa'QP ! P 

P - - l + «P /5 

Pj^-Z-Pq is the quantity Eegnault has designated by l+100a, let us denote 
it by A ; then if, for convenience, we put 

r - l+100a', f = a+lOOa", 

each determined value of A will give the equation of condition 
y + P^.y' = A + AP P.a. 
The following are Regnault's determinations of A augmented, in general, 
by 0.00018, for the reason that we adopt the mean coefficient 0.00018153 for 
the expansion of mercury between 0° and 100°, found by this experimenter, 
instead of the value & ^ So used by him, (See Note, p. 31 of the vol.) ; the last 
column contains the page of the volume, where the experiments may be found. 

P . P x . A. No.Obs. Obs.-Cal. Page. 



Q m 


.110 


m 


.149 


1.36500 


10 


—0.00012 


99 





.174 





.237 


1.36531 


3 


—0.00004 


99 





.266 





.362 


1.36560 


2 


—0.00003 


99 





.375 





.510 


1.36598 


4 


+0.00005 


99 





.548 





.746 


1.36673 


3 


+0.00038 


57 





.756 





.7535 


1.36724 


4 


+0.00035 


66 





.557 





.754 


1.36651 


18 


+0.00014 


43 





.656 





.757 


1.36641 


14 


—0.00022 


33 





.747 


1 


.016 


1.36663 


3 


—0.00014 


58 





.771 


1 


.049 


1.36696 


11 


+0.00014 


51 


1 


.678 


2 


.286 


1.36778 


2 


—0.00059 


109 


1 


.693 


2 


.306 


1.36818 


4 


—0.00021 


109 


2 


.526 


2 


.517 


1.36962 


2 


—0.00018 


114 


2 


.622 


2 


.614 


1.36982 


2 


—0.00011 


114 


2 


.144 


2 


.924 


1.36912 


2 


+0.00007 


109 


3 


.656 


4 


.992 


1.37109 


4 


+0.00031 


109 
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Adopting, for convenience, as an unknown in the place of y, 
x= r +f— 1.367, 
we have the following equations, to each of which we attribute a weight 
equal to a tenth of the number of experiments it is founded upon : 

x — 0.7086/ — 0.3268a = — 0.00200 
x — 0.6064/ — 0.4398« = — 0.00169 
x — 0.4821/ — 0.5790a = — 0.00140 
x — 0.3534;-' — 0.7237a = — 0.00102 
x — 0.1728/ — 0.9258a = — 0.00027 
x — 0.1674/ — 1.1410a = + 0.00024 
x _ 0.1670/ — 0.9354a = — 0,00049 
x — 0.1650/ — 1.040 a = — 0.00059 
x + 0.0105/ — 1.131 a = — 0.00037 
x + 0.0317/ — 1.155 a = — 0.00004 
x + 0.708 / — 1.913 a = + 0.00078 
x + 0.718 / — 1.924 a= + 0.00118 
x + 0.818 / — 2.496 a = + 0.00262 
x + 0.863 / — 2.558 a = + 0.00282 
x + 1.004 / — 2.243 a = + 0.00212 
x + 1.834 / — 3.175 a = + 0.00409 
The derived normal equations, for determining x and /, are 
x— 0.0047/ — 1.162a = — 0.000144, 
— 0.0415a + 2.9547/ — 3.373a = + 0.007074, 
whence 

x =z — 0.000133+1.168a, 

/= + 0.002392+1.158a, 
and 

a! = + 0.00364475 +0.00010a, 

a"=^ 4- 0.00002392+0.00158a. 
The equation which determines a has already been obtained from the dis- 
cussion of the first class of experiments; it is 

f + 4 :l 4 l a '! = —0.0024286. 
1 + 4.747a' 

The last three equations being solved, we gather that the volume of any 
mass of air is represented by the formula 
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V= jri+aP^+(a'+a / 'P^)2'"], 



in which 

a = — 0.002565, a' = + 0.0036445, 

a" = 4. 0.00001987, /? = 0.64776. 
How well the second class of experiments is satisfied by this formula 

may be seen from the numbers in the column headed Obs.-Cal. 

If we have r = -9gg^ = 129°.l, V takes the form 

P 

Hence we have the noteworthy result that: — 

About the temperature 130°, air follows quite exactly the law of Mariotte. 
For the following temperatures and pressures the volume vanishes : — 

T. P. 



0° 


9995 m .49 


— 50 


4420 .13 


—100 


2048 .00 


—150 


896 .26 


—200 


314 .23 



These numbers may be regarded as indications of the magnitude of pres- 
sure necessary for the condensation of air. The table is in accordance with 
the well-known fact that reduction of temperature facilitates condensation. 

A table is given below which will be found useful in the application of 

the formula. It contains the functions log(l +aP") and -— — =-=-, the lat- 
ter being the coefficient of expansion under a constant pressure. 

As an example let us suppose that the volume of a mass of air has been ob- 
served under the pressure 2 m .5 and the temperature 20° ; it is required to 
find the factor necessary for reducing it to the pressure m .76 and tempera- 
ture 0°. From the table we get 3.07064. By employing the ordinary 
formula with the coefficient 0.003665 of expansion, there is obtained 3.06- 
482, which differs from the preceding by about a 525th part. 

Eigorously observations of pressure made in localities having an intensity 
of gravity different from that which prevails at Regnault's laboratory ought 
to be multiplied by the ratio of the former to the latter. The latitude of 
Eegnault's laboratory is stated at 48° 50' 14", the altitude above sea level 
at 60 m , and the intensity of gravity at 9 m .8096. 
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P. log(l +aPP). Coeff.ofExp. 



P. log(l+ «P/*). Coeff.ofExp. 



O m .O 
.1 
.2 
.3 
.4 
.5 
.6 
.7 
.8 

.9 

1 .0 
1 .1 



.2 
.3 



1 .4 

1 .5 

1 .6 

1 .7 

1 .8 

1 .9 

2 .0 

2 .5 

3 .0 

3 .5 

4 .0 

4 .5 

5 .0 

5 .5 

6 .0 
6 .5 



0.000000 
9.999749 
9.999607 
9.999489 
9.999384 
9.999288 
9.999199 
9.999115 
9.999035 
9.998958 
9.998885 
9.998814 
9.998745 
9.998678 
9.998613 
9.998549 
9.998487 
9.998426 
9.998367 
9.998308 
9.998251 
9.997979 
9.997725 
9.997485 
9.997257 
9.997039 
9.996829 
9.996626 
9.996430 
9.996239 



251 

142 

118 

105 

96 

89 

84 

80 

77 

73 

71 

69 

67 

65 

64 

62 

61 

59 

59 

57 

272 

254 

240 

228 

218 

210 

203 

196 

191 

186 



0. 



0036445 
36511 
36548 
36579 
36607 
36632 
36656 
36678 
36698 
36719 
36738 
36756 
36775 
36793 
36809 
36826 
36842 
36858 
36874 
36889 
36904 
36977 
37044 
37106 
37166 
37224 
37280 
37333 
37385 
37436 



66 
37 
31 
28 
25 
24 
22 
20 
21 
19 
18 
19 
18 
16 
17 
16 
16 
16 
15 
15 
73 
67 
62 
60 
58 
56 
53 
52 
51 
49 



7 m .O 

7 .5 

8 .0 

8 .5 

9 .0 
9 .5 

10 .0 

10 .5 

11 .0 

11 .5 

12 .0 

12 .5 

13 .0 

13 .5 

14 .0 

14 .5 

15 .0 

15 .5 

16 .0 

16 .5 

17 .0 

17 .5 

18 .0 

18 .5 

19 .0 

19 .5 

20 .0 

20 .5 

21 .0 
21 .5 



9.996053 
9.995872 
9.995695 
9.995521 
9.995352 
9.995185 
9.995021 
9.994860 
9.994702 
9.994546 
9.994393 
9.994242 
9.994093 
9.993946 
9.993800 
9.993657 
9.993515 
9.993374 
9.993236 
9.993098 
9.992962 
9.992828 
9.992695 
9.992563 
9.992432 
9.992303 
9.992174 
9.992047 
9.991921 
9.991795 



181 
177 
174 
169 
167 
164 
161 
158 
156 
153 
151 
149 
147 
146 
143 
142 
141 
138 
138 
136 
134 
133 
132 
131 
129 
129 
127 
126 
126 



0, 



0037485 
37533 
37580 
37626 
37671 
37715 
37758 
37801 
37843 
37885 
37925 
37965 
38005 
38044 
38083 
38121 
38159 
38196 
38233 
38269 
38306 
38342 
38377 
38412 
38447 
38482 
38516 
3g550 
38584 
38618 



48 
47 
46 
45 
44 
43 
43 
42 
42 
40 
40 
40 
39 
39 
38 
38 
37 
37 
36 
37 
36 
35 
35 
35 
35 
34 
34 
34 
34 



ON AD J USTMENT FORMULAS. 



EY E. L. DE E0EEST, M. A. 

(Continued from page 86.) 

Three separate tests of good adjustment have been proposed by the pres- 
ent writer, for which see Smithsonian Report of 1873, p. 333, and a more 
recent essay in pamphlet form j {Interpolation and Adjustment of Series, 
New Haven, 1876). The first of these tests is, that if we denote by v the 
residual error, found by subtracting an adjusted term from the correspoding 



